Quantum dynamics of the Neel vector in the antiferromagnetic molecular wheel CsFe 8 
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The inelastic neutron scattering (INS) spectrum is studied for the antiferromagnetic molecular 
wheel CsFe8, in the temperature range 2 - 60 K, and for transfer energies up 3.6 meV. A qualitative 
analysis shows that the observed peaks correspond to the transitions between the L-band states, 
from the ground state up to the S = 5 multiplet. For a quantitative analysis, the wheel is described 
by a microscopic spin Hamiltonian (SH), which includes the nearest-neighbor Heisenberg exchange 
interactions and uniaxial easy-axis single-ion anisotropy, characterized by the constants J and D, 
respectively. For a best-fit determination of J and D, the L band is modeled by an effective 
SH, and the effective SH concept extended such as to facilitate an accurate calculation of INS 
scattering intensities, overcoming difficulties with the dimension of the Hilbert space. The low- 
energy magnetism in CsFes is excellently described by the generic SH used. The two lowest states 
are characterized by a tunneling of the Neel vector, as found previously, while the higher-lying states 
are well described as rotational modes of the Neel vector. 

PACS numbers: 75.10. Jm,71.70.-d 



I. INTRODUCTION 

Molecular nanomagnets have recently attracted enor- 
mous attention, as these perfect magnetic nanoclusters 
were demonstrated to exhibit a number of spectacular 
quantum phenomena. In systems such as the single- 
molecule magnets Mni2 and Fes, quantum tunneling of 
the magnetization, and related effects such as quantum 
interference, were observed iS^i^ These phenomena are 
intimately connected to a high-spin ground state, so that 
at low temperatures these clusters essentially behave like 
large single spins S, where S = 10 for Mni2 and Fes- 
This single-spin (or giant-spin) model is in fact extremely 
successful in describing the magnetism of single-molecule 
magnets^ The actual underlying many-spin nature of the 
system, despite providing a large S, is not key to the ob- 
served tunneling phenomena. 

Another class of interesting molecules are the antifer- 
romagnetic (AFM) wheels, in which the magnetic cen- 
ters are arranged into almost perfect rings (only wheels 
with an even number of equivalent spin centers shall 
be considered) i 6 i 7 i 8 i 9 A variety of quantum phenomena 
were observed, such as magnetization steps ^i^i^i^i^ii^ 
quantized rotation of the Neel vector or quantized spin 
wavesji^ and quantum tunneling of the Neel vectorii§ii£ 
In AFM wheels, in contrary to single-molecule magnets, 
the ground state has a total spin of S — and hence is 
nonmagnetic. A single-spin model (with S = 0) would 
thus be inappropriate as it would not allow one to de- 
scribe the interesting spin dynamics in these systems, 
which arc manifestly many-spin effects. The lowest-lying 
excitations can be described in terms of the Neel vector, 
n, which is essentially a vector parallel to the magnetiza- 
tion of one of the AFM sublattices [in a classical picture 
it could be for example parallel to the up-pointing spins 
as in Fig. Hfl))]. 

The magnetism of the AFM wheels is very well de- 



scribed by the generic spin Hamiltonian 
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which includes isotropic (Heisenberg) nearest-neighbor 
exchange interactions, characterized by the coupling con- 
stant J < 0, and an uniaxial magnetic anisotropy, of the 
easy-axis type, described by an on-site anisotropy with 
D < 0. 18 N is the number of spin centers, Sj is the spin 
operator of the zth ion with spin s, and z is the uniaxial 
anisotropy axis perpendicular to the plane of the wheel. 

At small magnetic anisotropy, the lowest-lying excita- 
tions are characterized by a (quantized) rotation of the 
Neel vector, which in the quantum mechanical energy 
spectrum shows up as a low-lying rotational mode, the 
L band. This band is comprised of the lowest state in 
each S sector, i.e., the lowest S = 0, 1, 2, ...states, 
and their energies increase quadratically with S accord- 
ing to E(S) cx S(S + i),9 1 i9 1 2Q 1 2i,22 At hig^r energies, 
further rotational modes exist, collectively denoted as E 




n = +z 



FIG. 1: (Color online) (a) Crystal structure of CsFes [with 
Fe as dark gray spheres (red online) and H atoms omitted], 
(b) Classical ground-state spin configuration for the Neel vec- 
tor [long up arrow (dark green online]) oriented in the +z 
direction. 
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band, which correspond to the quantized spin-wave exci- 
tations expected in AFM spin systems Ai2*22i This picture 
of the excitations has been demonstrated in great detail 
experimentally by inelastic neutron scattering (INS) on 
the AFM wheel Cr 8 (for which J = -1.46 meV, D = 
-0.038 meV, or S /h = Ns^/TdJI = 2.7) 

If the magnetic anisotropy becomes strong, and ex- 
ceeds the threshold So/h > 4j2i the situation changes 
completely for the lowest lying excitation, which is 
then characterized by quantum tunneling of the Neel 
vector This scenario has recently been observed 
by INS in the title compound, the AFM wheel CsFe 8 
[the structure is shown in Fig.^a), here J = -1.80 meV, 
D = -0.05 meV, or S /H = 4.7] .±2 

In this work, the AFM wheel CsFes is further inves- 
tigated by INS measurements, extending the previous 
work, Ref. Qj| to a larger temperature and energy trans- 
fer range, which allows us to probe a much larger part 
of the low-lying energy spectrum. In particular, the L- 
band states have been detected from the ground state up 
to the S — 5 multiplet at an energy of ca. 15 meV. This 
provides a detailed insight into the dynamics of the Neel 
vector in the CsFes wheel. 

In the next section II, the experimental details and the 
data are presented. The analysis of the data has been 
split into two parts. Initially, a qualitative analysis is 
provided in section III. Subsequently, in section IV some 
preparatory work for a quantitative analysis is developed, 
with the final analysis presented in section V. Section VI 
provides a discussion of the results, and the manuscript 
concludes in section VII. 



II. EXPERIMENTAL 
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FIG. 2: Inelastic neutron scattering spectra of CsFeg at T 
= 2.2 K, 9.7 K, and 17 K, measured with A = 3.8 A. Panel 
(a) shows the neutron energy-loss, and panel (b) the neu- 
tron energy-gain side of the spectrum. The question mark 
indicates a weak feature which will be identified as a further 
transition VII' later on. 



The CsFes material, which has the chemical formula 
[CsFe 8 L 8 ]Cl • 2CHC1 3 • 0.5CH 2 C1 2 • 0.75L • 2.5H 2 0, 
with L = N(CH2CH20)3, was prepared by the same pro- 
cedure as in Ref. 1271 but crystallized from a 1:1 mix- 
ture of CHCI3 and CH2C1 2 by pentane vapor diffusion. 
The compound crystallizes in the monoclinic space group 
Pna21, and the molecules exhibit approximate C4 sym- 
metry with Fe-Fe distances of 3.142 - 3.164 A, Fig.GJa). 
The INS spectrum of CsFes was measured with the time- 
of-flight spectrometer IN5 at the Institute Laue-Langevin 
(ILL), Grenoble, France. 4 g of a non-deuterated pow- 
der sample was sealed under a helium atmosphere in 
a double- walled hollow Al cylinder (50 mm in height, 
16 mm external diameter, and 2 mm thickness) and in- 
serted in a standard ILL orange 4 He cryostat. INS spec- 
tra were recorded at temperatures between 2.2 K and 
60 K for incident neutron wavelengths of A = 3.8 A and 
5.0 A (energy resolution at the elastic peak was 161 [ieV 
and 121 /«eV, respectively). The data were corrected for 
detector efficiency using a vanadium standard; no further 
corrections were applied. The data shown correspond to 
the sum over all detector banks; the error bars are smaller 
than the symbol size. 



Figure [5] presents the neutron energy-loss and energy- 
gain sides of the INS spectrum as measured at A = 3.8 A 
for the temperatures T = 2.2 K, 9.7 K, and 17 K. FigureEl 
presents the A = 5.0 A data obtained for the tempera- 
tures T = 2.4 K, 9.7 K, 17 K, and 60 K. 



III. QUALITATIVE ANALYSIS 

The interpretation of the INS spectra is relatively 
straightforward. In the 3.8 A neutron energy-loss spec- 
trum at 2.2 K, Fig. Efa), two peaks I and II at about 
0.5 meV and 1.3 meV, respectively, are observed. They 
correspond to transitions from the ground state to two 
next higher-lying levels at the corresponding excitation 
energies, as indicated in Fig. Ufa). At 9.7 K the inten- 
sities of peaks I and II decrease, and a new transition 
III at 2.3 meV develops, which is assigned to a transi- 
tion from the first excited level to a higher-lying level at 
2.8 meV. Some additional weaker features at 1.5 meV 
and 2.0 meV (peaks IV and V) are visible, which grow 
stronger at 17 K. They are assigned to transitions from 
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FIG. 3: Inelastic neutron scattering spectra of CsFeg at T = 
2.4 K, 9.7 K, 17 K, and 60 K, measured with A = 5.0 A. Panel 
(a) shows the neutron energy-loss, and panel (b) the neutron 
energy-gain side of the spectrum. 



the second excited level to higher-lying levels at again 
2.8 meV and 3.3 meV. At 17 K a further strong feature 
at 2.9 meV, peak VI, is observed. As it is essentially 
absent in the data at the lower temperatures, it is as- 
signed to a transition from the 2.8 meV level to a level at 
about 5.7 meV. On the neutron energy-gain side of the 
data, Fig.JSfb), all corresponding anti-Stokes peaks (la- 
beled with an additional apostrophe) are observed with 
the expected temperature dependence. This confirms the 
above assignment of the levels, also demonstrating the 
magnetic origin of the transitions. 

The 5.0 A neutron scattering data shown in Fig. 
confirm the findings from the 3.8 A measurements with 
higher resolution (the 2.4 K and 9.7 K data were commu- 
nicated previously in Ref. Ilq) . The 60 K spectrum pro- 
vides valuable additional information. On the neutron 
energy-loss side, Fig. E^a), resolved transitions are no 
longer observed; however, an increased inhomogeneous 
intensity as compared to the 2.4 K data is evident. Ap- 
parently, at 60 K, a large number of transitions with 
widely varying transition energies become possible, and 
individual transitions cannot be resolved. We interpret 
this as the onset of transitions into the so-called quasi- 
continuum in the spectrum of AFM molecular wheelsAS 
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FIG. 4: (Color online) Energy spectrum of CsFeg as (a) de- 
rived from the neutron scattering data, and (b) calculated 
numerically from the microscopic Hamiltonian H, eq. 
using the best-fit parameters determined in this work. The 
level energies are given in meV. The labeling of the transi- 
tions, the classification of the levels, and the transitions a 
and P axe discussed in the text. 



In the neutron energy-gain spectrum at 60 K, Fig. EST)), 
the anti-Stokes peaks I' to VI' are clearly visible, as well 
as two further transitions VII' and VIII' at 3.9 meV and 
4.8 meV, respectively (the presence of peak VHP will be 
unambiguously confirmed in Fig. [B] below). Very weak 
indications of peak VII' can be found in both the 3.8 A 
and 5.0 A neutron energy-gain data at 17 K. This peak 
is thus assigned to a transition from the excited 5.7 meV 
level to a further higher-lying level at about 9.6 meV. 
Based on the weaker intensity of peak VIII', it is as- 
signed to a transition from the 9.6 meV level to a level at 
about 14.4 meV. The magnetic origin of peaks VII' and 
VIIF, and their assignment, is less evident from the data, 
but is confirmed by the quantitative analysis presented 
in section IV. 

With the above thermal-population arguments the en- 
ergy spectrum as compiled in Fig. Ufa) is directly con- 
cluded from the INS measurements. For the further anal- 
ysis, the following procedure was applied to each data 
curve on the neutron energy-loss side. The background 
was approximated by an analytical function and sub- 
stracted from the data. The corrected spectrum was then 
least-squares fitted to an appropriate number of Gaussian 
peaks. The decomposition into the different contribu- 
tions is explicitly shown for the 3.8 A data in Fig. [3 the 
fitting results for the peak positions, linewidths, and in- 
tensities are compiled in the Tables given in appendix A. 
The temperature variation of the peak intensities are in 
accord with the expected thermal Boltzmann population, 
supporting the above analysis. 

It is interesting to inspect the linewidths provided in 
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FIG. 5: Analysis of the A = 3.8 A neutron energy-loss spec- 
tra for T = 2.2 K, 9.7 K, and 17 K. The data curves for T 
= 9.7 K and 17 K were shifted for clarity by 0.1 and 0.175 
arb. units, respectively. Open circles represent the data and 
solid lines fits, which included an approximated background 
(dotted lines) and a number of Gaussian lines (dashed lines). 
The parameters of the Gaussian lines are given in Table A. II. 



Tables A.I and A. II. The linewidths of transitions I, II, 
and V are consistent with the experimental resolution 
for both the 3.8 A and 5.0 A data sets. This indicates 
that the energy levels at 0.5 meV, 1.3 meV, and 3.3 meV 
correspond to single states in the energy spectrum of 
CsFes. The linewidths of transitions III and IV, which 
are about 0.21 meV (FWHM) for both peaks, are con- 
sistently larger than the experimental resolution. This 
indicates that the level at 2.8 meV in Fig.^Ja) actually 
relates to two (or more) close-lying states in the energy 
spectrum. Finally, the linewidth of transition VI, ca. 
0.34 meV (FWHM), is even larger, which indicates that 
also the level at 5.7 meV contains two or more close- lying 
energy states. These findings are represented in Fig.^Ja) 
by the thicker gray bars. 

A similar analysis was performed for the 5.0 A neu- 
tron energy-gain data at 60 K in an energy range which 
embraces peaks VI', VIF, and VIII'. The decomposition 
of the data is shown in Fig. HJJ the best-fit values for the 
three Gaussian lines are given in Table A. III. The in- 
set of Fig. confirms the statistical significance of peak 
VIII'. The observed linewidths of 0.33-0.4 meV are sig- 
nificantly larger than the experimental resolution, with 
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FIG. 6: Analysis of a part of the A = 5.0 A neutron energy- 
gain spectra at T — 60 K. Open circles represent the data and 
the solid line the best fit, which included an approximated 
background (dotted lines) and three Gaussian lines (dashed 
lines). The parameters of the Gaussian lines are given in 
Table A. III. The inset shows the data after substraction of 
the background (black squares) and the three fitted Gaussian 
lines (solid line). 



similar conclusions as for the 3.8 A data. 

A qualitative assignment of the observed transitions 
is achieved by comparing the level scheme in Fig. EJa) 
with the energy spectrum as calculated numerically from 
H, eq. In order to calculate the energy spectrum, 
one could use the J and D values inferred previously 
from high-field torque measurements. 13 Since these val- 
ues, however, are very close to the best-fit values to 
be determined below in section V, the energy spectrum 
obtained for the best-fit values, which is presented in 
Fig. Elk), will be used here. In Fig. Efb), the states are 
labeled by the quantum numbers S and M of the total- 
spin operator S = J2i S» (using the notation \S, M) for 
the levels). However, it should be stressed here that S 
is not a good quantum number for CsFes, because of the 
rather large magnetic anisotropy present in the system^ 
For the moment, S should be taken as an arbitrary addi- 
tional label for distinguishing the energy levels; a further 
detailed discussion of the situation will be given in sec- 
tion VI. 

From a comparison of Figs.^a) and^fb), the identi- 
fication of the experimentally observed transitions is ap- 
parent, considering the INS selection rules AS — 0, ±1 
and AM = 0,±1. The 0.5 meV and 1.3 meV levels 
correspond to transitions from the ground state |0, 0) 
to the levels 1 1,0) and |1,±1), respectively. As demon- 
strated in Ref. JJj, peak I reflects the quantum tunnel- 
ing of the Neel vector characterizing the low-energy spin 
dynamics in CsFes- It was also noted that the transi- 
tion between the first and second excited levels (denoted 
as a) is not observed, although it is formally allowed 
by the INS selection rules. The intensity of a is very 
weak compared to that of the other transitions, since the 
sublattice-magnetization vectors are of mesoscopic size 
in CsFes. Peaks III, IV, and V correspond to transitions 
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from the 5 = 1 states to those of the 5 = 2 multiplet. 
Interestingly, peaks III and IV are confirmed to consist of 
two close transitions, in agreement with the observed en- 
larged linewidths (the transition |2, ±1) — * |2, ±2), which 
is expected at ca. 0.6 meV at higher temperatures, turns 
out to be too weak in intensity to be detectable) . Peak VI 
corresponds to transitions from the 5 = 2 states into the 
5 = 3 multiplet. The zero- field-splitting of these multi- 
plets, which allows six close transitions, is reflected by the 
broad linewidth of peak VI. Finally, peaks VII and VIII 
correspond to transitions from the 5 — 3 to the 5 = 4, 
and the 5 — 4 to the 5 = 5 multiplcts, respectively. 
Here again, the spread in the transition energies of the 
individual transitions due to the zero-field-splitting is too 
small to be resolved experimentally, hence the broadened 
peaks. 

The energy level diagram shown in Fig. EJa) is one of 
the most extended level schemes determined for a molec- 
ular nanomagnet from INS measurements, yet all essen- 
tial features of the data can be readily explained by basic 
arguments. The above analysis, in our opinion, provides 
a nice textbook example for the analysis of INS data from 
molecular clusters. 



IV. EFFECTIVE HAMILTONIAN AND 
INELASTIC NEUTRON SCATTERING 
INTENSITY 

In the next section a more sophisticated analysis of 
the data will be presented, which yields best-fit values 
for J and D, as well as some insight into the accuracy 
and limitations of the generic spin Hamiltonian A 
fit of the data to INS spectra calculated from the mi- 
croscopic Hamiltonian using standard least-squares 
fitting routines, however, is numerically expensive for a 
system as large as CsFes. Another procedure has thus 
been developed, the essentials of which are outlined in 
this section. 

In our experiments, only transitions between states 
of the so-called L band were observedi2ii£i2£ The L 
band, also known as the tower of states^ or the quasi- 
degenerate joint statesjSSiSi arises from a quantized rota- 
tion of the Neel vector in AFM wheels and other systems. 
It consists of the energetically lowest-lying states for each 
5, whose energies increase as E(S) cx 5(5 + 1). This 
set of states can be very well described for small AFM 
Heisenberg rings by approximating the wave functions 
by IPaPbSaSbSM), with S A = S B = 10 for CsFe 8 (S A 
and Sb denote the total spin of sublattices A and B, re- 
spectively. Pa and [3b abbreviate intermediate quantum 
numbers, but are omitted in the following) Physi- 
cally, this approach works well because the internal spin 
structure due to the dominant Heisenberg interaction is 
essentially classical J''- 5 In this approximation, the effec- 
tive (two-sublattice) Hamiltonian 



with J = 0.5366J and D = 0.1870D is obtained for 
CsFes ^ Since Hab is designed to describe only the 
states of the L band, numerical calculations become much 
cheaper (for CsFes , the dimension of the Hilbcrt space is 
1679616 while that of H AB is only 441). 

Concerning the energies, Hab has been demonstrated 
to yield very accurate results^ The calculation of INS 
spectra, however, needs further consideration since not 
only the energies but also the local transition matrix ele- 
ments and so-called interference factors need to be prop- 
erly approximated. While transition matrix elements are 
expected to be reproduced reasonably well, the situa- 
tion is less obvious for the interference factors. In fact, 
these factors, which typically produce a pronounced os- 
cillatory behavior of the INS scattering intensity as func- 
tion of momentum transfer (Q dependence), arise from a 
correlation of the transition matrix elements and the ge- 
ometrical structure of the molecule^S^i Hab essentially 
replaces the octanuclear spin-5/2 ring structure of CsFes 
by a dimer of two spin- 10 centers, which is obviously a 
very different geometrical structure. A naive calculation 
of the INS spectra from Hab should thus be suspected 
to produce incorrect Q dependencies, and hence incorrect 
INS intensities. For CsFes this is indeed the case. In the 
following, however, it will be shown that the interference 
effects can be retained by a suitable generalization of the 
INS scattering formula. The procedure is actually very 
simple. 

For a spin cluster, the INS cross section is given by^S 
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(3) 

where C{Q,T) = (je 2 /m e c 2 )(k' /k) exp[-2IF(Q, T)] 
(where all symbols have the usual meaning), = 
l/(fcsT), and Z(T) is the partition function. For a pow- 
der sample in zero magnetic field one finds^i 

I nm (Q) = J^F*(Q)F j (Q){^j (QR ij )S i -S j 

+i 2 (Q%)E T i 2) *(^) T i 2) (S^i)},(4) 



where Fi(Q) is the magnetic form factor of the ith 
spin center, j)~ is the spherical Bessel function of or- 
der k, and Ry = R^ — Rj is the distance vector be- 
tween the zth and jth ion. An explicit expression for 

Eq. is given in Ref. |3 Tq (v) is the gth component 
of the spherical tensor of rank k constructed from the 
Cartesian vector v, and T g (S^Sj) represents the ten- 
sor product [TW(Si) <g> T( 1 )(S i )]^ 2) . The ordered prod- 
ucts Si a Sjp, which appear in the explicit expression of 

Tq 2 \SiSj), stand for (n\Si a \m)(rn\Sjp\n), where \n) de- 
notes an eigenstate of the Hamiltonian under considera- 
tion. Equation (0| can be also written more compactly 
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WQ) =^2^2fif(Q,Kij)U^\S i S j ), (5) 

ij kq 

whereby introducing the interference factors f^ q (Q, Ry) 
(k = 0, 2 and \q\ < fc) and the symmetrized spherical 
tensors Uq^ (which are proportional to Re[Tg fc ^] for q > 

( k) 

and Im[Tg ] for q < 0, see appendix B). 

The easiest way to derive an effective formula for the 
INS intensity is to follow the procedure suggested in 
Ref. |3f| for the construction of an effective Hamiltonian 
for AFM spin clusters with essentially classical internal 
spin structure. It starts by noting that a classical spin 
structure relates to a mean-field situation for each of the 
AFM sublattices. Accordingly, all spins Sj on a sublat- 
tice v are replaced by the mean-field spin S u /N u , where 
S„ = ^2iev ^ ano - N v is the number of spins on sublat- 
tice v. 

Concerning the INS intensity, one replaces all matrix 
elements Sj a with i £ v by S ua /N v , puts all S ua outside 
the brackets, and immediately obtains 



Vfl kq 



with the effective interference factors 
1 



fui(Q) 



N U N, 



■ E 



(6) 



(7) 



In the case of CsFes, with its two sublattices A and B, 
Eq. |[BJ yields a scattering formula analogous to that of 
a dimer, but with the dimer interference factors fn(Q) 
and /i2(Q) replaced by the effective interference factors 
Jaa{Q) and Jab{Q)- In the scattering formula, now only 
the matrix elements Sau and Sbu appear, which can be 
calculated from the eigenpairs obtained from Hab- 

Figure |7| presents numerical results for the Q depen- 
dence of the INS intensity in CsFes as calculated from 
the microscopic spin Hamiltonian H (solid lines), the ef- 
fective spin Hamiltonian Hab (dotted lines), and the ef- 
fective Hamiltonian Hab but now using Eq. (|HJ) for the 
INS intensity (dashed lines) [in the first two cases Eq. (QJ 
was used] . Noting Fig. [7| it is clear that the pure dimer 
model is not capable of producing correct Q dependen- 
cies. This is not surprising, and just reflects the fact that 
a dimer model can only describe the interference between 
two ions, while in the full cluster interference occurs be- 
tween eight ions. This also explains why the exact Q 
dependence shows more oscillations per Q range. In con- 
trast, the above procedure of reinstating the full geomet- 
rical structure of the molecule works very well. For the 
intra-multiplet transition a (|1,0) — * |1,±1)), the exact 
and effective dimer curves are essentially indistinguish- 
able in Fig. [7{c). The curves in Figs. |7fa) and (b) also 
coincide with each other if one scales the effective dimer 
curves by factors of 0.940 and 0.911, respectively. These 
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FIG. 7: Q dependence of the neutron scattering intensity of 
CsFe8 as calculated from the full Hamiltonian £Q| (solid lines) , 
the dimer Hamiltonian Hab (dotted lines), and the dimer 
Hamiltonian Hab using the effective INS formula Eq. © 
(dashed lines) . Results are shown for the three transitions (a) 
|0,0) [1,0), (b) |0,0} -► |1,±1), and (c) [1,0) |1,±1>, 
see Fig. [1] In the calculations with the pure dimer model, the 
distance between the two spins Sa and Sb was set to 3.2 A, 
corresponding to the iron- iron distance in CsFe8- 



numbers also give an idea of the accuracy of the transi- 
tion matrix elements obtained from the effective Hamil- 
tonian approach, which apparently is better than 10%. 
Furthermore, Fig. [7| clearly shows that for Q > 0.5 A the 
intensity of the intra-spin-multiplet transition a is indeed 
two orders of magnitude weaker than that of transitions 
I and II as noted previously^ 



V. QUANTITATIVE ANALYSIS 

In this section a complete quantitative analysis of 
the experimental data is provided. Typically, spin- 
Hamiltonian parameters are obtained by least-square fit- 
ting the calculated energies (or transition energies) to the 
peak positions obtained from a Gaussian analysis of the 
data, similar to that shown in Fig. However, since 
in the present case some of the observed peaks are com- 
posed of several individual but close-lying transitions the 
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experimental spectra were fitted directly by simulated 
INS spectra. 

Only the results of the analysis for the 3.8 A neu- 
tron energy-loss data at 17 K will be shown explicitly, 
as this data set contains the most information and thus 
allows the best judgment as to the agreement between 
experiment and theory. A similar analysis for the other 
data yields redundant results. For each pair of values for 
J and D, theoretical INS spectra were calculated from 
the eigenpairs of Hab using the effective INS formula 
Eq. For the final best- fit values, and some other 

cases, spectra were calculated from the microscopic spin 
Hamiltonian H, using Eq. in order to confirm the 
accuracy of the results. In the calculations, the crystallo- 
graphic positions of the iron(III) ions were used, and the 
form factors Fi(Q) were estimated by the standard an- 
alytical approximations^ The variation of the Q range 
with energy transfer hio of the IN5 spectrometer was ac- 
counted for by numerically integrating I(uj,Q), where 
Q = Q(uj,9), over the scattering angle spanned by the 
detectors. Q(u>,6) is related to the energy transfer frui 
and the angle 9 of the scattered neutron beam via the 
scattering triangle. 

The result of a least-squares fit to the 3.8 A, 17 K neu- 
tron energy-loss data is shown in Fig. Eta) • The agree- 
ment is excellent, with the calculated curve reproduc- 
ing all the details present in the data. The obtained 
best-fit parameters were J — -1.80(2) meV and D = 
-0.050(1) meV, which are in agreement with previous 
findings The linewidth was set as a free parameter 
in the fit, yielding 0.18 meV, very close to the nominal 
instrumental resolution of 0.161 meV. The good agree- 
ment between experimental and calculated curves is also 
notable because several transitions contribute to the var- 
ious peaks, as indicated by the vertical solid lines in 
Fig. |Et a )- Peaks I, II, and V correspond to individual 
transitions, while peaks III and IV are each made up of 
two close-lying transitions, and peak VI is the result of 
six close-lying transitions, in perfect agreement with the 
conclusions drawn from the analysis of the linewidths in 
section III [see also Fig. Ufa)]- The explicit assignment 
of all the transitions which form the spectrum in Fig. [S] 
is given in Ref. |37j From the expected transitions be- 
tween the S = 0, 1, 2, and 3 multiplets, all but two 
are observed. These are the two intra-multiplet tran- 
sitions |2,±1) -> |2,±2) and |1,0) -> |1,±1) expected 
at 0.571 meV and 0.818 meV, respectively. Their inten- 
sity is two orders of magnitude weaker as compared to 
the other transitions, because, as mentioned already and 
discussed in Ref. Ilo . the sublattice-magnetization vectors 
are of mesoscopic size in CsFes . The other intra-multiplet 
transitions, while allowed, occur at too low energy and 
are thus buried under the elastic line. 

Using the best-fit parameters, the transitions between 
the S = 3, 4, and 5 multiplets can also be calculated. 
The result agrees nicely with the peak positions of peaks 
VII and VIII, thus confirming their magnetic origin, and 
their assignment as transitions from the S = 3 to S = 4 
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FIG. 8: (a) Best-fit curve, to the background-substracted 
data, for the A = 3.8 A neutron energy-loss spectrum for T 
= 17 K (J = -1.80 meV, D = -0.050 meV). The vertical lines 
indicate the transitions and their intensities which contribute 
to the neutron spectrum, (b) Calculated neutron energy-loss 
spectra. The solid curve was calculated with the best-fit pa- 
rameters using the effective Hamiltonian approach, Eq. ©, 
the short-dashed curve was calculated with the best-fit pa- 
rameters using the full Hamiltonian, Eq. @, and the dashed 
curve was calculated from the full Hamiltonian with a 10% 
distribution in J and D (see text for details). 



and S = 4 to S = 5 multiplets, respectively. 

As a summary of the analysis so far, the transitions 
between the S = 0, 1, 2, 3, 4, and 5 multiplets (of the L 
band) were observed and could be reproduced very well 
by the Hamiltonian H. Apparently, the generic Hamil- 
tonian H, which contains only the two free parameters 
J and D, provides an excellent description of the excita- 
tions in CsFes, within the experimental resolution of the 
present experiments of about 0.1 meV. 

Indications have been found in previous works on other 
AFM wheels that Hamiltonian (JJJ, though being an ex- 
cellent approximation, is not complete and needs to be 
extended by further weak termsiii*i^ii22iiSi22i This raises 
the question as to how well Hamiltonian (JTJ in fact de- 
scribes CsFes. It should be noted that for some wheels, 
effects which are due to coupling of the spin system to 
the environment have been found. 40,41 however these arc 
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a different subject and not of interest here. The pres- 
ence of further weak terms seems particularly relevant 
for CsFes as it does not exhibit eight-fold symmetry; the 
molecule investigated in this work has a crystallographic 
C2 axis. Accordingly, two deviations from the perfect sit- 
uation described by Hamiltonian JJJ were considered, i) 
an additional biaxial anisotropy term and ii) a variation 
of the exchange and anisotropy constants along the ring. 

Biaxial anisotropy contributions for rather symmetric 
wheels were recently deduced for the Cr§ and the Feio 
wheelsiiiSi In these cases, the microscopic Hamiltonian 
Hx = H + £Ei(3i,x - S? y ) was used. The Cr 8 wheel 
exhibits a crystallographic C4 symmetry (though indica- 
tions of a structural phase transition at low temperatures 
seem to have been found), however the INS data were 
significantly better described with \E\ = 0.004(1) meV, 
corresponding to a ratio \E/D\ = 0.11^4 For the Feio 
wheel, which exhibits crystallographic C2 symmetry, \E\ 
= 0.0036(2) meV was found, corresponding to a quite 
substantial ratio \E/D\ = 0.21ill 

In the case of CsFes, an upper limit for the value of | ^E7| 
can be inferred from the present data as follows. A bi- 
axial anisotropy would split the two |1, ±1) levels, which 
are degenerate in a uniaxial model. Accordingly peak 
II should split into two [see also Fig. Ufa)]. The 5.0 A 
data at low temperatures yields a linewidth for peak II of 
ca. 0.13 meV (see Table A. II), very close to the nominal 
experimental resolution. Thus, a splitting of, for exam- 
ple, 0.08 meV (corresponding to an estimated linewidth 
of peak II of ca. 0.18 meV) would have been clearly de- 
tected in this experiment. Since the splitting of the |1, 0) 
and |1, ±1) levels by the D term is ca. 0.8 meV, one thus 
can conclude that \E/D\ < 0.1, or with the determined 
value for D that \E\ < 0.005 meV. For a more quantita- 
tive analysis, the 3.8 A, 17 K neutron energy-loss data 
was least-squares fitted to the effectivc-Hamiltonian ver- 
sion of .Hi, with J, D, and E (and the linewidth) as free 
parameters. Independent of the starting values, the fit 
converged to the same values for J and D as above, and 
to \E\ < 5 X 10~ 5 meV, that is, values which are zero 
within the statistical error. The inclusion of a biaxial 
term thus did not improve the quality of the fit. 

A more suitable technique for detecting a possible 
E term, or other additional terms in the Hamiltonian, 
is electron paramagnetic resonance (EPR), as it allows 
one to detect the transitions within a spin multiplct 
with higher resolution (in the case of CsFeg this is ca. 
0.01 meV compared to ca. 0.1 meV for the INS exper- 
iment). EPR experiments at 35 GHz and 190 GHz on 
single-crystals of CsFes were thus performed (and will 
be reported elsewhere^). The EPR spectra are well 
described using only uniaxial anisotropy, which sets a 
rather stringent upper limit on the strength of a biaxial 
anisotropy. Thus from a magnetic perspective CsFes is a 
very uniaxial and symmetric wheel. 

As a second possible deviation from the generic situ- 
ation, a variation of the exchange and anisotropy con- 
stant within the ring was considered, as described by the 



Hamiltonian H 2 = J i$i • S l+ i + J2i (there 

enforcing Y^i^i = J an d Si -A = D). A variation 
might originate, for instance, from the reduced C2 sym- 
metry of the ring, which gives rise to four slightly differ- 
ent Fe-Fe distances, or from disorder in the crystallites. 
In order to test such an effect, Ji and Di values were 
drawn from a Gaussian distribution with a 10% stan- 
dard deviation, and the INS spectra calculated. This is 
shown in Fig. |HJb) , and is compared to the calculation 
for the generic Hamiltonian (JJJ [Fig. Efb) also presents 
the INS spectrum as calculated with the effective Hamil- 
tonian approach, just to reinforce the high accuracy of 
this procedure]. It is apparent that a 10% variation of 
the parameters is quite possible within the present exper- 
iments. This is not unexpected because, as discussed in 
Ref. a change of the coupling and anisotropy param- 
eters along the ring has no effect on the energies of the 
states of the L band up to first order. However, this is 
not so for the states of the E band, which correspond to 
the quantized spin- wave excitations on the ring^SS Ob- 
servation of the spin waves in CsFes, expected at around 
7.5 meV and 10 meV at low temperatures, could thus 
allow further insight into this issue. 



VI. DISCUSSION 

In the previous sections it has been shown that the 
generic Hamiltonian H, Eq. ft}, which involves only the 
two parameters J and D, provides an excellent descrip- 
tion of a significant portion of the low-lying excitations 
in CsFes- In the analysis, the states were labeled by the 
spin quantum number S, which, however, since D is non- 
zero, is not an exact quantum number (though M still 
is). The question thus arises as to how good a quantum 
number S is. The situation for CsFes is in fact quite 
ambivalent. 

In our previous INS studyi^ it has been demonstrated 
that the lowest energy transition |0,0) — ► |1,0), or peak 
I, corresponds to the Neel-vector tunneling transition. 
Here, the influence of the magnetic anisotropy is so strong 
that S completely loses its significance as a quantum 
number for the two involved states. The label S only 
can be used in the sense that the states are connected 
adiabatically to the corresponding states at zero mag- 
netic anisotropy. This is visible in Fig. [HI which presents 
the calculated energy spectrum as function of the ratio 
D/J (using J = -1.80 meV). For small D/J the states 
obviously can be labeled by S and M ("small" means 
here D/J < 0.02, Ref. S3). For large D/J, where S cer- 
tainly is no longer a good quantum number, every state 
nevertheless can be uniquely related back to one of the 
\S, M) states at small D/J. Hence S can still be used as 
a label, although the actual wave functions by no means 

are reasonably approximated by eigenfunctions of S . 

The loss of the significance of S with increasing D/J 
can be inferred from the energy spectrum, Fig.|3 at least 
qualitatively (for a more rigorous discussion one would 
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M=0,±2 



0.35 



FIG. 9: Low-lying energy spectrum of an octanuclear s = 5/2 
wheel as function of D/J, as calculated from H with J = - 
1.80 meV. At small anisotropics, the levels can be classified by 
the spin quantum numbers S and M. For large anisotropies, 
S loses its significance, though M remains an exact quantum 
number. The dashed line indicates the parameter D/J as 
obtained for CsFes. 



have to inspect the wave functions, but here a qualita- 
tive consideration shall suffice). For D/J = 0, the lev- 
els |5, ±M) are degenerate since 5 is an exact quantum 
number. For small values of D/J the levels split, but 
the states corresponding to a particular 5 are still clus- 
tered together. For large D/J, however, the pattern of 
the level clustering changes completely. For instance, at 
small D/J the gap between the 5 = ground state and 
the two S = 1 levels is much larger than the splitting 
between the two 5 = 1 levels. This indicates that 5 is 
a good quantum number and that the three states |0,0), 
|1, 0), and |1, ±1) are approximately eigenfunctions of 5 2 . 
For larger D/J, however, the two states |0,0) and |1,0) 
come close to each other, while |1, ±1) splits off and joins 
the |2, ±1) level. In this situation the wave functions are 

clearly no longer eigenfunctions of S . 

It is interesting to note in Fig.|5|that the value of D/J 
for this crossover apparently increases with the energies 
of the spin states. Concerning the states of the 5 = 
and 5 = 1 multiplets, the crossover is at about D/J — 
fj.021 2 ^ At this value, however, the states of the next- 
higher lying 5 = 2 multiplet (and the multiplets with 
5 > 2) are still clustered together; only at significantly 
higher D/J values are these levels split so strongly that 
they start to join levels coming from other spin multiplets 
(the |2, ±1) level approaches the |1, ±1) level of the 5 = 
1 multiplet, and |2,0) and |2, ±2) merge with levels of 
the 5 = 3 multiplet). From Fig. [5] it can thus be seen 
that with increasing magnetic anisotropy D/J the levels 
successively lose their 5 character, with larger crossover- 
D / J values for larger 5. 

This provides a rather intuitive picture of the Neel- 
vector dynamics in AFM wheels. The situation at zero 
magnetic anisotropy has been well discussed in previous 
works&i^iiSiSfliSL^ The lowest-lying excitations, which 
form the L band, may be characterized as a quantized 




FIG. 10: (Color online) Qualitative sketch of the semi- 
classical picture for the dynamics of the Neel vector in AFM 
wheels, (a) Small magnetic anisotropy, where the dynamics 
is characterized by rotation of the Neel vector, (b) magnetic 
anisotropy exceeding So /h > 4, where the low-energy dynam- 
ics is characterized by a tunneling of the Neel vector, and (c) 
larger magnetic anisotropy, where successively further higher- 
lying states are driven into the tunneling regime. 



rotation of the Neel vector (that is, in a classical picture 
the Neel vector may precess freely). 

The introduction of an easy-axis anisotropy means 
that (in a classical picture) the orientations of the Neel 
vector along the anisotropy axis (z axis) are energeti- 
cally favored over the other orientations, which is usually 
sketched by a double-well potential with two minima for 
the Neel vector along the z axis, and an energy barrier 
in between (see Fig. I1U0 . For small anisotropies the en- 
ergy barrier is small and the precession of the Neel vec- 
tor is largely undisturbed [Fig. ITUT a)] , and thus may be 
observed as for instance in the molecule Cr 8 ii£ With in- 
creasing anisotropy the energy barrier increases, so that 
the (classical) ground state eventually falls below the top 
of the barrier, and the dynamics is characterized by tun- 
neling of the Neel vector, Fig. [TUr b)i 24 i 25 i 26 In the quan- 
tum mechanical picture, this affects the 5 = and 5 = 1 
multiplets, which lose their 5 character. For the higher- 
lying multiplets, the anisotropy is not yet strong enough 
to compete with the higher rotational energies involved 
in these states (as indicated by their larger 5 values), 
so that they are still well characterized by a rotation 
of the Neel vector, and by 5. With further increasing 
anisotropy, the two lowest-lying states are driven deeper 
into the tunneling regime (as indicated in Fig. by the 
exponentially shrinking gap between |0, 0) and |1, 0)), but 
eventually also the next-higher levels become tunneling 
states (as indicated in Fig.[5|by the exponentially shrink- 
ing gap between |1,±1) and |2,±1)), see Fig. ^fc). This 
continues with further increasing anisotropy. 

Apparently for CsFes, the situation portrayed in 
Fig. llOf b') is realized, where the low-energy dynamics is 
characterized by a tunneling of the Neel vector (and the 
quantum number 5 has lost its meaning), while at the 
higher energies the dynamics can still be rationalized as 
a rotation of the Neel vector (so that the spin multiplets 
with 5 > 2 retain their 5 character). The 5 = 2 multiplet 
nevertheless shows an interesting signature of the strong 
anisotropy, namely the |2,0) and |2,±1) levels are very 
close to each other in CsFes (see Fig. [!|| ■ This indicates a 
strong fourth-order term B\ for the 5 = 2 multiplet with 
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a ratio -B4/D2 = 5 x 10 4 (for comparison in Mni2>^ 
B°/D 2 = 5x 10- 5 )ii 

VII. CONCLUSIONS 

In this work the previous INS measurements 1 — on the 
AFM wheel CsFes were extended to a larger temperature 
and energy transfer range. This has allowed us to observe 
the L-band states from the ground state up to the S = 
5 multiplet. The analysis of the data, on a qualitative 
level, was straight forward and of textbook quality. For 
the quantitative analysis the effective Hamiltonian con- 
cept, which is well-known to be very successful in the 
interpretation of energy spectra, was extended such as to 
also facilitate an accurate calculation of the INS scatter- 
ing intensities. This allowed a least-squares fitting of the 
data. The results showed that the low-energy magnetism 
in the CsFes wheel is excellently described by the generic 
spin Hamiltonian H , Eq. ft|. with only two magnetic pa- 
rameters. The role of the quantum number S has also 
been discussed qualitatively, and showed that in CsFes 
the situation is ambivalent: On the one hand, the very- 
lowest excitations are characterized by a tunneling of the 
Neel vector, implying that S is no longer a good quantum 
number, while on the other hand the higher-lying states 
are still well described as rotational modes of the Neel 
vector, for which S is a good quantum number. 

It would be clearly interesting to synthesize and iden- 
tify AFM wheels with an anisotropy so strong that the 
situation of Fig. llOf c) is realized, i.e., where also higher- 
lying states are driven into a tunneling regime. The semi- 
classical theory for such a situation has not been, to the 
best of our knowledge, yet developed. It would be, how- 
ever, important to put the qualitative picture of the Neel- 
vector dynamics in strongly anisotropic wheels outlined 
in section VI on a sound theoretical basis. 

For very large D/J, the situation should bear a close 
analogy to that found in the (Mn^ molecule, which rep- 
resents a dimer of two Mn4 single-molecule magnets cou- 
pled by a weak AFM interaction^Si^Li^i This molecule 
attracted quite some interest recently because of its dis- 
tinct magnetic hysteresis curve. According to the effec- 
tive two-sublattice Hamiltonian Hab, Eq. (J2J, an AFM 
wheel can be considered magnetically as a dimer of two 
large spins Sa and Sb- Thus, for small J, and because 
of the easy-axis anisotropics of the spins Sa and Sb , one 
formally would have exactly the situation of a dimer of 
single-molecule magnets with a weak AFM interaction, 
similarly as in (Mn^. It is thus predicted that in highly 
anisotropic AFM wheels phenomena such as exchange- 
biased magnetic hysteresis curves could be observable. 
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APPENDIX A 

As described in section III, the experimental INS spec- 
tra were analyzed by approximating the background with 
an analytical function and least-square fitting of the cor- 
rected data to an appropriate number of Gaussian lines. 
The resulting parameters for the six peaks I to VI vis- 
ible in the 3.8 A energy-loss spectra [Figs. Ufa) and |5] 
are listed in Table A.I, those for the peaks I, II, and IV 
visible in the 5.0 A energy-loss spectra [Fig. Efa)] are 
listed in Table A. II, and those for the three peaks VI' 
to VHP visible in the 3.8 A energy-gain spectra at T — 
60 K [Figs. Elb) and © are listed in Table A.III. The 
results of a similar analysis for the other neutron energy- 
gain data, while in agreement, is not given here because 
of redundancy. The linewidths reported in the tables 
are the widths of the best-fit Gaussian lines (no convo- 
lution with instrumental resolution). For comparison, 
also the theoretically expected instrumental resolutions 
of IN5 are given in the last columns. These, however, are 
not achieved in the experiment (probably due to sample 
inhomogeneity). This is evident from the experimental 
resolution at the elastic line, which for A = 5.0 A (3.8 A) 
was 121 /xeV (161 ^eV) as compared to the theoretical 
instrumental resolution of 105 fieV (141 /zeV). Hence, for 



TABLE I: Results of the analysis of the A = 3.8 A neu- 
tron energy-loss spectra for the temperatures 2.2 K, 9.7 K, 
and 17 K. Each peak in the data was least-square fitted by 
a Gaussian line; the peak position, linewidth (FWHM), and 
intensity (in arb. units) are listed below. For further details 
of the analysis see section III and Fig. |S| The last column 
reports the theoretically estimated instrumental resolution of 
IN5 at the respective transfer energies. 



peak 


T 


peak energy 


width 


intensity 


instr. width 






(meV) 


(meV) 


(xlO -2 a.u.) 


(meV) 


I 


2.2 K 


0.5274(3) 


0.1732(8) 


1.31(1) 


0.132 




9.7 K 


0.5275(4) 


0.1747(9) 


0.761(4) 


0.132 




17 K 


0.5237(8) 


0.186(2) 


0.393(3) 


0.132 


II 


2.2 K 


1.3450(3) 


0.1647(6) 


1.382(5) 


0.119 




9.7 K 


1.324(3) 


0.141(5) 


0.5(1) 


0.119 




17 K 


1.327(4) 


0.140(7) 


0.24(5) 


0.119 


III 


2.2 K 


2.289(9) 


0.17(2) 


0.029(3) 


0.105 




9.7 K 


2.263(2) 


0.210(4) 


0.41(1) 


0.106 




17 K 


2.253(4) 


0.221(6) 


0.33(1) 


0.106 


IV 


9.7 K 


1.43(2) 


0.20(2) 


0.4(1) 


0.118 




17 K 


1.46(2) 


0.22(1) 


0.37(5) 


0.117 


V 


9.7 K 


2.036(4) 


0.134(6) 


0.09(1) 


0.109 




17 K 


2.042(3) 


0.140(4) 


0.13(1) 


0.109 


VI 


17 K 


2.937(2) 


0.338(3) 


0.354 


0.097 
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TABLE II: Results of the analysis of the A = 5.0 A neu- 
tron energy-loss spectra for the temperatures 2.4 K, 9.7 K, 
and 17 K. Each peak in the data was least-square fitted by 
a Gaussian line; the peak position, linewidth (FWHM), and 
intensity (in arb. units) are listed below. For further details 
of the analysis see section III. The last column reports the 
theoretically estimated instrumental resolution of IN5 at the 
respective transfer energies. 



peak 


T 


peak energy 


width 


intensity 


instr. wic 






(meV) 


(meV) 


(xlO -2 a.u.) 


(meV) 


I 


2.4 K 


0.5178(5) 


0.133(1) 


2.53(2) 


0.094 




9.7 K 


0.5211(2) 


0.1255(4) 


1.361(5) 


0.094 




17 K 


0.5203(5) 


0.12(1) 


0.66(1) 


0.094 


II 


2.4 K 


1.3171(4) 


0.1303(8) 


2.41(1) 


0.080 




9.7 K 


1.3040(5) 


0.110(2) 


0.84(4) 


0.080 




17 K 


1.303(1) 


0.101(4) 


0.30(3) 


0.080 


IV 


9.7 K 


1.398(5) 


0.184(4) 


0.86(4) 


0.079 




17 K 


1.420(5) 


0.208(5) 


0.80(3) 


0.078 



TABLE III: Results of the analysis of the A = 5.0 A neu- 
tron energy-gain spectra at T = 60 K. Each peak in the data 
was least-square fitted by a Gaussian line. The peak position, 
linewidth (FWHM), and intensity (in arb. units) are listed 
below. For further details of the analysis see sections III and 
IV, and Fig. Oil The last column reports the theoretically esti- 
mated instrumental resolution of IN5 at the respective trans- 
fer energies. 



peak 


peak energy 


width 


intensity 


instr. width 




(meV) 


(meV) 


(xlO -2 a.u.) 


(meV) 


VI' 


-2.945(4) 


0.363(8) 


0.112(3) 


0.187 


VII' 


-3.883(8) 


0.400(2) 


0.062(3) 


0.218 


vnr 


-4.83(4) 


0.33(9) 


0.010(3) 


0.252 



the neutron energy-loss side the experimental resolution 
should be assumed to be essentially independent on en- 
ergy transfer and of the value at the elastic line. The 
errors in the tables are the statistical errors returned by 
the fitting routine. For a complete error analysis instru- 
mental uncertainties (widths of time channels, inaccurate 
sample position, etc.), which amount to a few fJ-eV, also 
have to be considered. 

For a given peak, the best-fit values for the transition 
energies obtained from the various curves show some 
variation, however, this is well within the typical limits. 
For transition II, the obtained linewidths at the higher 



temperatures, see Tables A.I and A. II, are unrealistically 
small, however, this is not significant in view of the close 
proximity of transitions II and IV. 



APPENDIX B 

In this appendix, explicit expressions for the sym- 
metrized spherical tensors Uq k ^ and the interference fac- 
tors f kq are given, which were introduced via Eq. © in 
section IV in order to provide a compact notation for the 
INS scattering intensity of powder samples. The sym- 
metrized spherical tensors read 

Uq ^(SiSj) = Si X Sj X + SiySjy + Si Z Sjz, 

Uq (SjSj) = (2Si Z Sj Z Si X Sjx SiySjy)/\ // 6, 

U2 (SiSj) = (SixSjx SiySjy)/y/2 } 
U_ 2 (S^Sj) = SixSjy ~\~ SiySjx, 

u[ \SiSj) = SixSj Z + Si Z Sjx, 

U_i (SjSj) = SiySjz + SizSjy, (Bl) 

and the corresponding interference factors 

2 



fij {Q> ^1, . 



f!f 2 (Q,^ij) 
fij (Q, Rij) 



FtFj-joiQRij), 
F*F 3 j 2 (QRv) 31 ' n: 



R 



FtFjhiQRa 
F*Fjj 2 (QRij 
FfFjhiQRij) 



V6RI 

d2 _ r>2 

R 2 



FtFMQR^^jS 1 ^ 



(B2) 



In the previous literature, the form factors Fi(Q) were 
not included in the definition of the interference factors, 
which is physically sound, as these two factors have differ- 
ent physical origins. Here, however, they were included 
in the terms, because the approach described in sec- 
tion IV for the calculation of the INS intensity, from an 
effective Hamiltonian, is also very useful in the case of 
heteronuclear systems (such as the single-molecule mag- 
net Mn4Br, Ref . l48|) . In these cases the equations sim- 
plify considerably in the chosen notation. 
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